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Abstract 

GAIA will provide observations of some multiple asteroid 
and dwarf systems. These observations are a way to de- 
termine and improve the quantification of dynamical pa- 
rameters, such as the masses and the gravity fields, in 
these multiple systems. Here we investigate this problem 
in the cases of Pluto's and Eugenia's system. We simu- 
late observations reproducing an approximate planning of 
the GAIA observations for both systems, as well as the 
New Horizons observations of Pluto. We have developed 
a numerical model reproducing the specific behavior of 
multiple asteroid system around the Sun and fit it to the 
simulated observations using least-square method, giving 
the uncertainties on the fitted parameters. We found that 
GAIA will improve significantly the precision of Pluto's 
and Charon's mass, as well as Petit Prince's orbital ele- 
ments and Eugenia's polar oblateness. 

1 Introduction 

Astrometric monitoring of multiple systems is a powerful 
way in the Solar System to have access to the physical 
properties of small bodies. Indeed, the satellite motions 
provide the mass of the primary, as well as the harmonics 
of its gravity field. The bodies involved can be very differ- 
ent in size and masses, from dwarf planets to small aster- 
oids. These systems can be very compact and, as a result, 
are difficult to observe from Earth without adaptive optics. 
We can expect GAIA to observe the components of such 
systems, both the primary and its satellites, if these later 
are far enough from the primary Banceli n et al., 2011J . 
GAIA will make precise and regular observations of them, 
and as a result, will probably improve our knowledge of 
their dynamical parameters. 

The purpose of this paper is to estimate the precision 
on the dynamical parameters we can expect for multiple 
systems thanks to GAIA. This kind of systems have been 
discovered among nearly every family of small bodies. We 
have investigated the contribution of GAIA's observations 
for one system in the Kuiper Belt : Pluto, and one in 
the Main Belt : 45 Eugenia. After a presentation of the 
dynamical model used especially to describe multiple sys- 
tems in Section 1, we will develop the case of each system 



in a different section. 



2 Dynamical model 

We use here the same numerical model which has been de- 
veloped in |Beauvalet et al., 2012| . We consider the mo- 
tion of every bodies of a multiple system in the inertial ref- 
erence frame ICRF centered on the barycenter of the Solar 
System. We compute the motion of the bodies disturbed 
by the Sun and the planets, whose positions are obtained 
through the numerical ephemeris DE405 [Standisli7l 998 . 
The initial positions, velocities and masses in Pluto's sys- 
tem come from |Tholen et al., 2008] . When needed, we 
include the second order harmonics of the polar oblate- 
ness of the primary, J^. We use the following notations 



• i an integrated body from the considered system, 

• j the Sun or a planet, 

• mi the mass of the body i, 

• Yj the position vector of the body j with respect to 
the Solar System barycenter, 

• fjj the distance between bodies i and j, 

• Ri the equatorial radius of body I, J% the polar 
oblateness of body I, 

• Ujt the potential of the I body's oblateness on the i 
body's center of mass. 

We then obtain the following equation of motion : 

in 



Guij (r j — rj 



Grm(Yi - n) 



(1) 



+GmiViU n - GmiViU, 



where Ujt is a function of the oblateness of I, and 4>i the 
latitude of i with respect to I equator : 



l > 2 



(2) 



Numerical integration of the equations of motion has been 
made using the 15th order Gauss-Radau integrator devel- 
oped by [Ev erhart, 1985| . Then we can adjust our model 
to the observations through a least-square procedure with- 
out constraints. 
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Table 1: Characteristics of Pluto and its satellites 

3 Pluto's sytem 
3.1 Description 

In 2013, Pluto will be about 33 AU from the Sun. We 
consider a four bodies system : Pluto, its most massive 

Nix and 



satellite Charon |Christy and Harrington, 1978 



3.2 Data simulation 

Our goal is to determine the precision on the masses we 
can expect thanks to GAIA's observations. To do so, we 
simulate data at the moment of already existing and ex- 
pected future observations of the system. We then fit our 
model to the simulations and extract the 1-a uncertainty 
from the least-square method. We did not include noise 
to our simulation. This comes from the fact that we only 
want the statistical uncertainty, a quantity which depends 
only on the uncertainty of the observations, the influence 
of the parameter on the system and the correlations be- 
tween the parameters. We used the Rendez-vous software 
to obtain a possible schedule of observations of Pluto's 
sytem by GAIA between 2013 and 2017. This software 
code has been developed by C. Ordcnovic, F. Mignard 
and P. Tanga (OCA) for Gaia DPAC. The uncertainty at- 
tached to these simulated GAIA observations is considered 
to be 1 mas. As a first approximation, we neglect the fact 
that this precision is available only in the direction of the 



Table 2: Dates used for the simulation of GAIA observa- 
tions for Pluto's system. 

Dates 



Hydra jWeaver et al., 2005 . Physical details on these ob- 
jects are given in Table [T] A last satellite has been dis- 
covered in July 2011 Showalt er et al., 201 1] and has not 
been included in the model. 

The particularity of Pluto among other systems in the 
Kuiper Belt is that it will be observed in situ by the 
probe New Horizons in 2015. This means that the sys- 
tem will then be observed simultaneously by this probe 
and by GAIA. Nonetheless, New Horizons will observe the 
four bodies of the system during a very short amount of 
time, whereas GAIA will only be able to detect Pluto and 
Charon, but its observations will be regularly made during 
five years. The contribution of New Horizons' observations 
to our knowledge of the dynamical parameters of the sys- 
tem has been investigated in [Bcauvalc t et al., 2 012 . In 
this previous study, we have found that the only parame- 
ters which can be estimated are the masses of the bodies, 
while the oblate gravity fields will not be obtained, even 
at the time of New Horizons arrival. We use here the same scan 
method for GAIA. 
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The dates used for GAIA simulations are given in Ta- 
ble [2j The simulated ground-based observations consist 
of ten observations per year, with the same uncertain- 
ties as current ground-based observations. The New Hori- 
zons simulations are obtained using a preliminary sched- 
ule of the mission, as well as its estimated uncertainties. 
More details on these two sets of simulations are given in 



Beauvalet et al., 2012 



As a result, we have two different sets of simulated data 



1992-2014+NH : reproducing the existing observa- 
tions, the future possible observations and New Hori- 
zons temporary observation schedule 



• GAIA 
GAIA 



reproducing observations of the system by 



Simulations with GIBIS suggest that the two bodies, Pluto 
and Charon, should always be detected as separated ob- 
jects. Nix and Hydra have respective magnitude of about 
23.7 and 23 |Stern et al., 2006 , so these two bodies will 
not be detected by GAIA. 
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3.3 Results 



3.3.1 Contribution of GAIA's whole mission 



The uncertainties on the masses for every object of the 
system is given in Table [31 We can see that GAIA's ob- 
servations will lower the uncertainties of the masses of 
Pluto and Charon when comparing the cases with and 
without GAIA. We can also see that the precision of the 
masses of Nix and Hydra will also be a little bit modified. 
This comes from the fact that, if a parameter is more con- 
strained, the fitting process can no longer reduce the resid- 
uals through changing this parameter as much as before. 
As a result, the role, and so the uncertainty, of the other 
parameter will also change. The masses which will be the 
most improved by GAIA are Pluto's and Charon's. We 
could have expected that the uncertainties on the masses 
of Nix and Hydra would have been much lowered because 
of the stronger constraints put on Pluto and Charon. This 
is not what happens here, only the uncertainty on Hydra's 
mass being lowered. This comes from the fact that the cor- 
relations between the parameters do not decrease linearly 
with the number of observations available. 

3.3.2 Orbit enhancement before New Horizons 
arrival 

GAIA will be launched two years before New Horizons' 
arrival, and so will collect data before the fly-by We have 
searched whether GAIA would enable us to put stronger 
constraints on the masses and, as a result, on the body 
ephemerides. We have used the same method as before 
with simulations spanning between 2013 and 2015 before 
New Horizons' fly-by The set of simulations named 1992- 
2015 uses the dates of the currently available observations, 
the future observations before New Horizons' arrival and 
the observations from GAIA before 2015. 

The obtained uncertainties on the masses and the satel- 
lites' semi-major axis are given respectively in Tableland 

SI 

As can be seen, Charon's mass would be considerably 
improved by those observations. Concerning Nix and Hy- 
dra's dynamics, the most interesting results concern their 
semi-major axis, whose uncertainty would be lowered even 
though they are not observed by GAIA. We can deduce 
that these observations, if they were to be available in time 
to prepare New Horizons' arrival, would be very precious 
to constrain not only Pluto's and Charon's motion, but 
those of Nix and Hydra as well. We can see here that 
in this case, strengthening the constraints on Pluto and 
Charon also constrains Nix and Hydra. 



4 Eugenia's system 
4.1 Description 

(45) Eugenia is one of the few known triple asteroids. The 
primary, Eugenia, is far from spherical, its shape being ob- 
viously oblate from high resolution observations and light- 
curve inversion method. Its two satellites, Petit Prince and 
S2004(45) are quite close to their primary and far smaller 
than it. The outermost satellite, Petit Prince, is the first 
discovered satellite of the system Merline et ah, 1 999 
in 1999, while the second one, closer to Eugenia, has 
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Table 3: 1-er error bars on the masses given by least 
square method using different sets of simulated obser 
vations, using rap = 870.3 
km 3 s -2 , m N = 0.039 km 3 s 
s- 2 |Tholen et ah, 2008| . 



km 3 s 2 , m c = 101.4 
2 and m H = 0.021 km 3 



set of simulated observa- 1-er error bars on the semi- 
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major axis (km) 




Charon Nix Hydra 
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5.8 23 155 


1992-2015 


3.25 10 43 



Table 4: 1-a error bars on the semi-major axis using two 
different sets of simulations. The current estimation of the 
semi-major axis is ac — 19570.45 km, = 49242 km and 
a H = 65082 km |Tholen et ah, 2008| 
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Table 5: Characteristics of Eugenia and its satellites 



semi major axis (km) 


Petit Prince 


S2004(45) 
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Table 6: 1-cr error-bar on some of Eugenia's and its satel- 
lites' dynamical parameters using different sets of simula- 
tions 



been discovered in 2007 |Marchis et al, 2007| . Petit 
Prince semi-major axis is only 3% of Eugenia's Hill ra- 
dius, meaning both satellites are deep inside Eugenia's 
gravitational well. Yet, most recent studies of the sys- 
tem |Marchis et al., 2008[ |Marchis et al, 2010| imply that 
they have a non-negligible inclination with respect to Eu- 
genia's equatorial plane. From Eugenia's shape, a theo- 
retical value of its second order polar oblateness has been 
estimated to be 0.19. This is in contradiction with a much 
lower value of 0.06 deduced from the satellite orbital mo- 
tions. S2004(45) is always too close to Eugenia to be ob- 
served by GAIA considering the difference in their magni- 
tude, about 8.2. On the contrary, Petit-Prince should be 
detectable when being close to its higher separation from 
Eugenia but the difference in their magnitude, about 7.4, 
prevents Petit Prince to be seen by GAIA when too close 
to Eugenia. 

4.2 Data simulation 

We used again the simulation of asteroids'transit on 
GAIA's CCD with Rendez-vous software, and extracted 
the scheduled dates for Eugenia's observations. The dates 
used are given in Table [71 As Petit-Prince will not be 
seen when too close to Eugenia, we decided to reject the 
positions where he would be closer than 500 mas from Eu- 
genia. We gave the observations a 1 mas uncertainty, as 
for Pluto's system. 

4.3 Results 

The uncertainties obtained for the semi-major axis of the 
satellites, and the pole orientation of Eugenia, are given 
in Table [B] Concerning the second order polar oblateness, 
its uncertainty for the 1998-2010 set is 0.0006, and 0.0002 
with GAIA observations, to be compared with J2 = 0.060. 
As can be seen when comparing the two sets of observa- 
tions, the orbital elements of Petit Prince will be con- 
strained by GAIA observations, and as a result, the un- 
certainties on those of S2004(45) will also be lowered. We 
also find that the pole orientation will be more constrained 
since this data is obtained from the satellites'motion. 



5 Conclusion 



Table 7: Dates used for the simulation of GAIA observa- 
tions for Eugenia's system. 
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For both systems, Pluto's and Eugenia's, GAIA will put 
new constrains on the dynamical properties of the system's 



4 



bodies. Pluto has the advantage that it will be observed 
by both GAIA and New Horizons, and each mission will 
give us new highlight on the system. The most interesting 
feature is that even though GAIA will not observe every 
body of the systems, the properties of the non-detected 
objects will also be better known than before the mission. 
The work presented here concerns only two multiple sys- 
tems. Though the fainter components of each of them will 
not necessarily be detected, because of its absolute mag- 
nitude or relative magnitude to the primary and/or its 
closeness to the primary, GAIA will give us new informa- 
tions about their dynamical behavior, masses and gravity 
fields, and hence new informations about the composition 
of the system. 
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